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Abstract 

We study the effects of coupling a cosmologically rolling scalar field to higher 
order curvature terms. We show that when the strong coupling scale of the 
theory is on the 10- 3 -10" 1 ey range, the model passes all experimental bounds 
on the existence of fifth forces even if the field has a mass of the order of the 
Hubble scale in vacuum and non-suppressed couplings to SM fields. The reason 
is that the coupling to certain curvature invariant acts as an effective mass that 
grows in regions of large curvature. This prevents the field from rolling down 
its potential near sources and makes its effects on fifth-force search experiments 
performed in the laboratory to be observable only at the sub-mm scale. We 
obtain the static spherically symmetric solutions of the theory and show that 
a long-range force appears but it is turned on only below a fixed Newtonian 
acceleration scale of the order of the Hubble constant. We comment on the 
possibility of using this feature of the model to alleviate the CDM small scale 
crisis and on its possible relation to MOND. 
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If we want to obtain a scalar field model of dynamical Dark Energy (DE) its mass 
has to be of the order of the Hubble scalqj (Hq ~ 10~ 33 eV) so that it evolves in the 
appropriate time-scale. This, of course, raises the question of how to avoid conflicts 
with the stringent limits on the (non) existence of extra long range forces coming 
from laboratory experiments and Solar System observations. Even if the scalar field 
couples to ordinary matter only with non-renormalizable interactions suppressed by the 
Planck scale, laboratory limits imply that its mass should be larger than 10 _3 eV pQ. 
Several mechanisms such as approximate global symmetries to suppress its couplings 
[2J or an environment-dependent effective potential and mass (the so-called chameleon 
mechanism) [3 4J have been proposed to reconcile a cosmologically rolling scalar field 
with the negative results of fifth-force search experiments. In this letter we consider the 
effect of coupling the quintessence field to higher order curvature invariants. We will 
see that such couplings provide a new way to suppress the effects of the force mediated 
by the scalar field at the Solar System or laboratory level and, at the same time, 
offer a range of unique observational signatures detectable in laboratory experiments 
or astrophysical observations. Specifically, the coupling we will be interested in is 

where is the ultra-light field, R^xa ^ e Riemann curvature tensor and A some mass 
scale. Its origin could be in a coupling of the scalar field with the Gauss-Bonnet 
combination, like in the DE models of [5], but since the scalar curvature or Ricci 
tensor part of such a coupling would not play any role in what follows we omit them 
for simplicity. And differently from the models in [5] we will not use this kind of 
coupling to drive the acceleration. We will assume that the field <fi has an appropriate 
quintessence-like potential to drive the acceleration and this coupling can be simply 
interpreted as a curvature-dependent mass for the scalar field. Its effect will be to 
keep the field firmly locked at <ft = in regions of large curvature, e.g. near sources. 
However, on a cosmological scale the field will have a small mass at late times and 
will be able to roll down its potential. But for this to happen we must require that 
the field mass obtained from this operator is small enough today (Am, < Hq). On a 
cosmological scale the operator above will translate into a mass of order 

Am2 ~ § (2) 



1 We will use natural units throughout the paper so that c and K are set to one. 
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so the scale A has to satisfy A > Hq ~ 10~ 33 eV for the field to be able to roll cosmo- 
logically. On the other hand if we want the theory to pass the laboratory tests even if 
the field has non- suppressed couplings to matter we need its mass to be m s > 10~ 3 eV 
locally. In the neighborhood of a massive object like the Earth the Kretchmann scalar 
reads R tlll \ a R ,J,vX ' T = 48-^njP-, where G is Newton's constant and M is the mass of the 
object, so that near sources the field gets a mass 

Plugging in this expression the Earth's mass and radius and imposing that this mass is 
larger than 10~ 3 eV we get the bouncj^l A < 10~ 31 eV. So we see that these constraints 
leave open just a small window for the scale A of two orders of magnitude. And notice 
that this window is in an extremely small mass range of the order of the Hubble scale, 
10~ 33 — lCT 31 e^, so one might wonder whether this theory could make sense as an 
effective quantum field theory at all. But to obtain the correct strong coupling scale 
of the theory we have to take into account the proper (canonical) normalization of 
the fluctuations of the metric. It is easy to see that the coupling (pQ), in a weak field 
expansion in flat space in powers of the fluctuations of the metric will give rise to 
couplings like 

T* R ^ R = WmI + ° \ A 2 M 3 ) (4) 

where we have suppressed the indices, and h^ c ' represents the canonically normalized 
metric fluctuations: = rj^ + h$/M p . We see then that the strong coupling scale 
in this theory is A s ~ (AMp) 1 ^ 2 that will be in the range 1CT 3 — 10 _1 eV. So our 
scalar field will be strongly coupled to gravity, with a strong coupling scale close to 
the vacuum energy scalqj. And notice also that the mass of the field locally in the 
Earth is pushed close to, or might even get above, the cut-off scale of the theory. So, 
depending on the chosen parameters, we might get out of the range of validity of the 
theory when computing its implications for laboratory searches of fifth forces. But 



2 Over the distance scales associated with laboratory experiments we can take the field mass as 
constant. And as long as the density of our probe is less than the Earth's, the Earth's gravitational 
field will give the dominant contribution to the scalar mass. 

3 In fact it has been argued in [3] that taking a cut-off in our effective gravitational theory at the 
vacuum energy scale might be a way to understand the smallness of the cosmological constant (see 
however [7]). 



2 



in any case we can expect that, like in the chameleon models of [3], the mass of the 
field will decrease and the range of the force mediated by this field will increase if an 
experiment is performed in space offering a very characteristic experimental signature. 

But even if the field is massive close to sources it is nearly massless in regions of 
small curvature. It is therefore natural to ask what is the effect of the coupling (TjQ) for 
long range forces. In order to gain some intuition about the expected behavior of the 
spherically symmetric solutions we can present the following argument. We expect that 
generically a massive boson will mediate a fifth force appreciable when the distance to 
the source is smaller than its inverse mass. In these situations the effect of the mass 
term in the equation of motion is negligible as compared with the terms coming from 
the Laplacian (i.e. momentum scale). We can approximate then the equation to solve 
by V 2 = 0" + 2(p' /r = 0, yielding the usual <fi oc 1/r profile characteristic massless 
fields. But for longer distances the mass term is not negligible anymore with respect 
to the derivative terms in the equations and this suppresses exponentially the solution. 
This well-known Yukawa exponential suppression (<f) oc e~ msr /r) is the reason why 
massive particles mediate only short range forces, significant only when r < m~ l . But 
now our effective mass term reads Am 2 ~ 48 (73^) 2 in the Schwarzschild geometry, 
so it depends on the distance to the source. And according to the previous argument, 
if the coupling above gives the dominant contribution to the mass we can expect an 
appreciable effect only whenever 

Ar 3 

r < m s (r)- 1 ~ — . (5) 

Since now the mass grows at short distance this condition means that the force would 
only appear at long distances, r > r c , where r c , defined as 

GM 



A, (6) 



corresponds to a fixed gravitational accelerations a c = A. We will confirm this result by 
explicitly obtaining the static spherically symmetric solutions to the equation of motion 
for the field with this r-dependent mass term. The fact that the force mediated by 
this field appears only below a fixed Newtonian acceleration scale makes this model 



4 This dependence of the mass with the distance yielding a fixed Newtonian acceleration as the 
threshold for the appearance of new effects is the same one finds for the scalar mode of the metric in 



the purely gravitational model of [8]. 
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interesting from a phenomenological point of view, since it can produce the emergence 
of a critical acceleration in dark matter (DM) halos helping with the problems of 
conventional cold dark matter (CDM) models at small scales. 

Let us present here a brief review of these problems, that will serve as an extra 
motivation for the kind of phenomenology that the present model can offer. The reader 
uninterested in this discussion can skip directly to the next paragraph. It can be said 
that the comparison of the structure of galaxies and DM halos formed in numerical 
simulations (within the ACDM paradigm) with their real observational counterparts 
is far from satisfactory. The main discrepancies come from the smallest scales: the 
simulations yield too much DM density in the central part of collapsed structures and 
too many sub-structures within DM halos (see e.g. [9][T0] and references therein). This 
is the main cause of what is sometimes called the CDM "small scale crisis" , and it has 
prompted the consideration of other flavors of DM, such as self-interacting [10] , warm 
[Tf] or decaying DM [12]. However it is unclear if a DM self-interaction is compatible 
with the structure at other scales like cluster halos [13j. Also, the possibility of having 
warm DM is severely constrained by observations of Ly-a forest [14] and it does not 
seem compatible with the evolved structure seen at high redshifts, since it acts delaying 
structure-formation [15] (and the same constraints limit decaying DM scenarios [16J). 
But on top of this, there are certain systematic features of galactic dynamics that 
seem hard to explain in any conventional DM scenario. In particular the onset of 
the mass discrepancy or need for DM seems to be tightly correlated with a critical 
acceleration scalqj, such that when the Newtonian acceleration produced by any given 
source is above it, there is no evidence for DM and when it is below we need to 
assume the presence of unseen forms of matter to explain the dynamics [T7]. And this 
critical acceleration scale is very small, curiously of the order of the current Hubble 
constant. The above considerations, plus the cosmological constant problem, could be 
interpreted as hints of the existence of non-trivial physics within the "dark sector"^. 
And the mechanism presented here seems a good candidate to alter the physics of 
the infrared domain of gravity in the right direction. It can explain the appearance 
of a fixed acceleration scale in DM halos and has the potential to save some of the 



5 This empirical fact underlies the success of MOND, a phenomenological modification of Newton's 
potential for small accelerations that is able to predict with remarkable accuracy the rotation curve 
of spiral galaxies from its visible matter content [TB] . 

6 These motivations are of a theoretical nature for the DE, but purely observational for DM. 
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contenders of CDM like warm or decaying DM by modifying the physics of structure 
formation because it can be seen as an enhancement of gravity in low density regions. 

We shall continue now exploring the implications of this kind of interaction for long 
range forces. As usual we will assume a coupling of this field to ordinary matter such 
that in the non-relativistic limit it can be parametrized as 

AS mt = J (7) 

where p is the matter energy density and (3 a dimensionless parameter. The equation 
of motion for the scalar field reads then 

where we have defined Q as the Kretchmann scalar Q = R^ y x a R^ uXa . We are neglecting 
the potential for </> since we are assuming that we are in the regime where the coupling 
(CQ) gives the dominant contribution to the mass and the terms in this equation coming 
from the potential V(0) are negligible. To find a solution one has to specify Q(x), 
i.e the background geometry. In principle this should be done self-consistently, which 
means solving at the same time the coupled gravitational equations taking into account 
the back-reaction of the scalar field configuration in the background geometry. We will 
neglect now this backreaction and we will check later the conditions under which this is 
actually a good approximation. In this case the Kretchmann scalar outside a spherical 
mass distribution reads Q = 48 (^pr) 2 , so the equation that we have to solve in vacuum 
in order to find the static solutions is 

d 2 (j) 2d(j) ( GM\ 2 



+ 



dr 2 r dr 

The general solution to this equation takes the form 



48 U^J 0=O - (9) 



%-l/4(*) + 7=^/4W (10) 



where I n (z) and K n (z) are the modified Bessel functions of the first and second kind 
respectively. We have defined z = 2v ^ A/ and the C\ are the integration constants. 
Notice that z(r) is basically the ratio of the gravitational acceleration to the critical 
acceleration a c = A. To determine the integration constants we have to satisfy the 
appropriate boundary conditions at infinity and in the source. As we said the field will 
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go to zero as we approach the sourceQ- For large z (or small r) we can approximate 

I-i/a{z) oc e z / z l l 2 and Ky^z) oc e~ z /z 1 ^ 2 . So requiring that the solution is well 

behaved at small r (remember that z oc 1/r 2 ) we see that we are forced to take C\ = 0. 

Notice that these solutions are analogous to the familiar oc e ±msr /r solutions for 

massive fields. Indeed, the heuristic argument presented before is backed by these 

exact solutions since e~ z ~ e ~ ms ^ r , and °^ I-i/i{ z ) I 'V^ is the analogous to the 

unphysical solution oc e +msr /r for massive fields. But curiously the role of the 

boundary conditions at the source and at infinity is reversed now with respect to the 

conventional case of a field with constant mass. In the conventional case we take 

the field to be asymptotically in its vacuum, say at = 0, so the coefficient of the 

e +m s r j r goi^ion is taken to be zero. The coefficient of the e~ msT /r term is obtained by 

satisfying the boundary conditions at the source giving the source mass as the constant 

of proportionality (the solution would be = (3 6 " sr for a coupling to matter like 

the one in (J7])). Now the field is locked at = at short distances, forcing as we said 

C\ — 0, while we allow it to be rolling at infinity, where it will have an unspecified, 

time dependent v.e.v. that we will call 0oo. As r — > oo we get the limit 

/ A \ 1/4 / A \ 1/4 

tf 1/4W /^2.3->/ 8 r(5/4)(— ) ~1.58(— ) , (11) 

where T(x) is the Euler gamma function. So the statiqj solution taking into account 
this boundary condition at infinity will be 

v 0OQ3 1 / 8 (GM\ 1/A 1 ^ (GM\ l/ * 1 

As we said this solution is exponentially suppressed at small distances (for z > 1) but 
at large distances it asymptotes to 0^ with a 1/r profile as can be seen from the large 
distance expansion 

3 1 / 4 r(-l/4) fGM' 

+ r(i/4) { 

(13) 




7 Actually, inside the source the equilibrium value for the field will be displaced from zero because 
the source acts as a tadpole through the coupling ([7]), but the field will go to zero exponentially on a 
distance scale ~ mj 1 once we get out of the source. 

8 In reality the solution is not static because we are considering time-dependent asymptotic bound- 
ary conditions. But we are assuming here that the relaxation time of the system is much smaller than 
the time-scale associated with the change in boundary conditions (~ H^ 1 ), so this static solution is 
a good approximation at any given time. 
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Figure 1: The solution for the scalar field with r c = 1. The dashed line corresponds to 
the asymptotic 1/r dependence. 

And notice that although this solution depends on the source mass and asymptotic 
field value, it is independent of the coupling of the field to matter, eq.([7]). We also see 
how we recover the usual 1/r 2 falloff in the force caused by <fi when r ^> r c . In fig.l we 
present a plot of this solution and in fig. 2 a plot of its gradient, that would correspond 
to the force that a test particle would feel due to the profile of 4> i n this spherically 
symmetric solution. It is apparent how this force is switched on when z ~ 1 (or r ~ r c ) 
and gets exponentially suppressed at smaller distances. 

But we still have to check that the backreaction of the scalar field configuration 
on the background is actually negligible, so that taking Q = 48 (^r) 2 is a good 
approximation. Since the scalar field is exponentially suppressed for r < r c , we know 
that the backreaction is negligible in this regime and the solution is valid. We will check 
now under which conditions neglecting the backreaction for r > r c is justified. This 
can be quantified by looking at the ^-dependent terms in the gravitational equations. 
The relevant terms in the equations of motion for the metric read 



where Tffi is the usual energy-momentum tensor for <\>. We can now interpret the 
right hand side of these equations as an "effective energy-momentum density" , , 
that will have an influence on the geometry. So if the "energy-momentum" stored 
in this r.h.s. of the Einstein equations is small compared to M we can neglect its 
influence on the background, while if it was bigger it can have a significant impact on 




(14) 
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Figure 2: The gradient of the scalar field (i.e. the force) with r c = 1. The dashed line 
corresponds to the asymptotic 1/r 2 dependence. 



the long-distance geometry. We are interested thus on the integration of this quantity 
as compared to M at any given radius. As mentioned before it is obvious that for r < r c 
we can neglect because is exponentially suppressed. It is easy to check that 

asymptotically in our solution, when r ^> r c and ~ 1/r, the r.h.s. of the gravitational 
equations has the dependence 

r ,m „ fegj£! (15) 

In order to get an estimation of the order of magnitude of the integration of we 
will approximate <f> by this long distance solution for r > r c and by zero for r < r c . We 
get then 

«<■»' = / ~ / ~ ^f- ~ & (16) 

J J A 3 r 7 A 3 r4 A v ; 

Comparing this with M we see that in order to be able to neglect the backreaction we 
must have 0^ < (AM) 1/2 . 

It is also interesting to compare the usual Newtonian gravitational force with the 
long distance force that will appear when r > r c as a result of this profile for the scalar 
field. Their ratio, for r 3> r c , is 



/ty(r)/M p 



(17) 



GM/r (AM) 1 / 2 ' 

Notice that, when (3 is of order one, the condition that we need for the backreaction of 
the scalar field on the geometry to be negligible is precisely that this ratio is small. So 
if P ~ 1 and the extra force due to the scalar becomes of the order of the gravitational 
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one we can not neglect the influence of the scalar field configuration on the background 
geometry and our solution is not valicj^. And since (AM) 1 / 2 <C M p for any reasonable 
M, this ratio will be very large unless (p^ <C M p as well, so we can only use the 
solutions we have presented when the v.e.v. of the field at infinity is small compared 
to the Planck mass. But although the precise form of the solution might be different 
when the backreaction is not negligible, we can still expect that the solution will be 
suppressed for r < r c (since the backreaction is negligible in this regime) and will 
recover a oc 1/r profile at large distances, r»r c (since the whole mass term for the 
scalar will be negligible in this regime and the field will behave as effectively massless). 
Nevertheless, it would be interesting to explore in more detail the behavior of the 
solutions when the backreaction is not negligible. In this case the solution for the metric 
will suffer an appreciable modification (with respect to the Schwarzschild solution), and 
this brings hopes that a model along these lines could be used to build MOND-like 
modification of gravity that bypasses the need for CDM. When the backreaction is 
negligible and the solution for the metric is not significantly altered we can not hope 
to build a realistic model that does not need DM because the effects of DM are also 
felt in gravitational lensing, and this effect can not be achieved without modifying the 
actual solution for the metric at long distances (see e.g. [19]). 

To conclude, we have explored the effects of coupling the quintessence field to higher 
order curvature invariants and in particular we have focused in a coupling that can be 
interpreted as a curvature-dependent mass for the scalar field. On a cosmological scale 
the field mass is small so it can roll down its potential acting as dynamical Dark 
Energy. But as we approach sources this mass grows, driving the field towards zero, 
and we have found that if this coupling is strong enough it can shield the Solar System 
from long range forces mediated by this scalar. It is also worth mentioning that this 
effect also implies that no time variation of fundamental constants associated with this 
scalar field would be observed in local experiments 10 !. However, in the places where 
the field is rolling, basically inter-stellar or inter-galactic space, a time dependence of 
fundamental constants could be observed. So this kind of model has also the potential 
to reconcile the claims of the detection of a time evolution of fundamental constants 



9 However, if (3 ^> 1 we can still have an extra force stronger than gravity appearing at long distances 
while the backreaction is small, so we can trust our solution. 

10 Or in the early Universe, since the field mass would be much bigger than the Hubble constant at 
early cosmological times and the field would be settled in its high-curvature minimum at <f> = 0. 
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derived from quasar absorption spectra with conflicting bounds obtained locally, like 
those coming from the natural Oklo reactor (see [20] for a review). Furthermore, the 
model also offers the opportunity to be tested in table-top laboratory experiments, and 
like the chameleon models of [3] boasts the characteristic signature of mediating a force 
whose range would be larger on space, where the local curvature is significantly lower. 
Needless to say, this theory has many more implications for astrophysics because at long 
distances an extra force will appear, becoming noticeable only when the gravitational 
acceleration produced by sources fall below a given scale. This feature makes the model 
interesting from a phenomenological point of view, since it appears to have the right 
properties to help with the problems of CDM at small scales. Interestingly, the strength 
of this force depends on the source's mass and asymptotic value of the field, but it is 
independent of the coupling of the field to matter. However a more detailed study of 
the effects of this kind of interaction in the formation and evolution of structures is 
beyond the scope of the present paper. 
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